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Abstract 

We discuss the crossover from the quantum critical, z = 1, to the quantum 

disordered regime in high-Tc materials in relation to the experimental data on 

the nuclear relaxation, bulk susceptibility, and inelastic neutron scattering. In 

our scenario, the spin excitations develop a gap A-^l/^ well above Tc, which 

is supplemented by the quasiparticle gap below Tc. The above experiments 

yield consistent estimates for the value of the spin gap, which increases as the 

correlation length decreases. 
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The underdoped cuprates YBa2Cu306.637 YBa2Cu408, and Lai.85Sro.i5Cu04, display 
strikingly different magnetic behavior from that measured for the fully doped system, 
YBa2Cu307. In YBa2Cu307, the uniform spin susceptibility, xo, is essentially tempera- 
ture independent, and the product ^^TiT of the nuclear spin-lattice relaxation time and the 
temperature has monotonic temperature dependence; in the underdoped compounds Xo is 
temperature dependent while the temperature dependence of ^^TiT is non-monotonic. Re- 
cent measurements of the spin-echo decay time, T2G for YBa2Cu306.63 [H, together with 
the earlier measurements for YBa2Cu306.9 0, provide additional valuable information on 
such spin gap phenomena in the underdoped systems. 

Alternative physical origins of the spin gap have been proposed recently by Millis and 
Monien 0] and Sokol and Pines 0, hereafter SP. Millis and Monien argued on the basis of 
bulk susceptibility measurements that the spin gaps in YBa2Cu306.63 and Lai 85Sro.i5Cu04 
possess a different physical origin, with the gap in YBa2Cu306.63 arising from a magnetic 
coupling between adjacent CUO2 layers, while for Lai.85Sro.i5Cu04 it is attributed to the spin 
density wave fluctuations of a metal. On the other hand, SP applied scaling arguments to the 
analysis of the NMR experiments to determine the magnetic phase diagram of the cuprates. 
They argued that at high temperatures these materials, and the closely related system, 
YBa2Cu408, are in the quantum critical (QC), z = 1, regime (see characterized 
by a temperature independent ratio ^^TiT/T2g, and a linear dependence on T of ®^TiT 
PJ^, while the common physical origin of the spin gap behavior these systems exhibit at 
lower temperatures is the suppression of the spectral weight for frequencies, < A ~ c/^, 
characteristic of the quantum disordered (QD) regime. In the quantum critical regime 
0J^, LJ ~ ~ T and z = 1, while in the quantum disordered regime ^ saturates is 
the antiferromagnetic correlation length, lj characteristic energy scale, and z dynamical 
exponent). In the present communication we show that the spin susceptibility and neutron 
scattering measurements on these underdoped cuprates are consistent with the SP magnetic 
phase diagram, reconcile the low frequency spin dynamics measured by NMR with the 
results of neutron scattering experiments, and give estimates of the spin gap. A, the spin 
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wave velocity, c, and ^. 

We consider first the Knight shift and bulk susceptibility measurements. As shown in 
Ref. U^, in the QC, z = 1, regime, xo is linear in T, in agreement with numerical calculations 
as well as experiment for La2Cu04. In the renormalized classical regime, xo is finite at 
T = 0; on the other hand, in the QD regime, the susceptibility rapidly (exponentially for 
the insulator) decreases due to the presence of the gap in the spin excitation spectrum, 
and tends to zero as T ^ [0. As may be seen in Fig.||, at high temperatures Xo(^) 
displays the expected linear in T behavior for La2Cu04, Lai.85Sro.i5Cu04, YBa2Cu3 06.63, 
and YBa2Cu408 with somewhat varying slope. 

The crossover from the QC to the QD regime is clearly visible in the experimental 
data shown in Fig.|l], apart from the insulator. As discussed in SP, the temperature T*, at 
which ^^TiT starts to deviate from its linear behavior (T* ~ 200K for YBa2Cu306.63 and 
YBa2Cu408, and T* ~ 125K for Lai.85Sro.i5Cu04), signals the beginning of this crossover. 
This downturn, when xo rapidly (although not exponentially in the doped case) decreases 
for T ^ A, reflects the suppression of the low frequency spin excitations. As may be seen 
in Fig.0, the downturn in Xo(^) begins at nearly the same temperature as that at which 
GSrp^rp ]^ggj^g deviate from its linear behavior. It follows that, as expected for the QD 
regime, the spin gap influences the low frequency dynamics both near q = and q= (vr, tt). 
The experimental finding that the ^^O spin lattice relaxation rate in YBa2Cu306.63 [ll3 and 
YBa2Cu408 follows the temperature dependence of Xo(T) [(^^TiT)^^ ~ xo(T)] is thus 
consistent with our scenario. If one assumes, with SP, that the effect of the quasiparticles on 
the QNLo" model parameters is somewhat less for Lai.85Sro.i5Cu04 than in YBa2Cu306.63, 
and ascribes this to a smaller hole density in the former material, and also, as Dupree 
et al. [|12| argue, that YBa2Cu408 possesses a hole density comparable to that found in 
YBa2Cu306.8, we see that the temperature at which the downturn occurs increases with 
hole concentration, going from ~ lOOK for Lai.85Sro.i5Cu04 to ~ 240K for YBa2Cu306.63 
and YBa2Cu408. Moreover, since the magnitude of Xo(T) in the QC and QD regime may 
be expected to increase with hole doping, this trend would be likewise evident in Fig.|l], once 
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one assigns a hole density in Lai.85Sro.i5Cu04 similar to that found for YBa2Cu306.55, and 
takes its orbital Knight shift to be in the range 0.02% to 0.076%. 

A close examination of the data shows that two successive downturns are visible in Fig.|l|, 
with only the second corresponding to the superconducting transition (Fig.|I]). This supports 
the view that in the superconducting state of QD materials, two separate gaps are present; 
a charge gap which reflects the pairing of quasiparticles in the superconducting state, and 
a spin gap A~l/^, which reflects the QD behavior modified, in the superconducting state, 
by the pairing correlations. 

We turn next to the evidence for a spin gap and a cross-over from the QC to the QD 
regime found in the neutron scattering experiments on the YBa2Cu306+2; system. Consider 
first YBa2Cu306.6, where the measurements of Ref. show that the q-integrated (i.e. 
local) dynamical response function, x'Li^) = / '^^qx"(ct, t^), is a universal function of u/T 
|r3| which behaves as x'L ~ high temperatures (i.e. small uj/T). This behavior 

is consistent with the QC scaling found for this compound in the analysis of the NMR 
experiments Moreover, the value of Xl{^) = Txl{^) inferred from experiment is in 



good agreement with the universal function calculated using the expansion 0. For low 
temperatures, the suppression of x'L YBa2Cu306.63, which is seen as an abrupt deviation 
from the universal behavior of x'L ^ function of uj/T, is an indication of the spin gap 
phenomenon. From an analysis of their experimental results, Tranquada et al. [0 conclude 
that the magnitude of the spin gap for this material is A ~ lOmeV. 

Since at T = 0, x'Li^) — ^ uj<A, an alternative approach is to associate the maximum 



of dx'L{^)/duj with the spin gap magnitude, A. Rossat-Mignod et al. |jT5[ find that for 
YBa2Cu3 06.69 this maximum occurs at 16 meV and we take this to be the value of A for 
this compound. At lower hole doping, in the range 0.4 < x < 0.5 just above the metal- 
insulator transition at x = 0.41, according to Ref. 0, YBa2Cu306+x is expected to be in the 
QC region for T > 50K. While no spin gap was seen in the non-superconducting YBaCu06.4, 



Rossat-Mignod et al. |TB| have reported the existence of a small gap in the low temperature 
superconducting region of YBaCu06.45 and YBaCu06.5i. 



At higher hole doping, in the heavily-doped YBa2Cu306+x with x ~ 1, a value of the spin 



gap ~ 26meV has been reported by Rossat-Mignod et al. using unpolarized neutrons [jT5 
Mook et al. [|I6| have used polarized neutrons to study YBa2Cu307 and have reported that a 
low-frequency suppression of the signal, ascribed to the spin gap formation, may be present 
at ~ 35meV; they find as well a magnetic peak at 41meV. Inelastic neutron experiments 
in the intermediate region of oxygen concentrations 6.75 and 6.9 have also been reported in 
Ref. [|14|, but the measurements have been performed only at low temperatures, where the 
superconducting gap, which complicates such an analysis and may affect the spin gap, is 
already developed. We show in Fig.^ the values of A which have been reported for various 
oxygen concentrations. 

Our conclusion that the spin gap increases as the doping increases is somewhat unex- 
pected in view of previous analyses but may be easily understood in terms of the QNLcr 



model, where A = c/^, provided ^ increases faster than c as the doping decreases. One can, 
moreover, combine the results of neutron scattering and NMR experiments to obtain an 
estimate of the doping dependence of c and ^. 

We consider first YBa2Cu306.63, for which we estimate A = 12meV by interpolating be- 
tween the experimental results for YBa2Cu306.6 and YBa2Cu306.69- This value is somewhat 
less than that found by applying scaling arguments to the experimental results for T2g(T), 
shown in Fig.^. Applying the approach of [0, one obtains that in the QC regime, T2g{T) 
is linear in T, while the ratio of this slope and its extrapolated intercept is the product of 
a universal number, C,, and A. On taking C, as obtained from the leading order of a 1/N 
expansion, one finds A = 18meV; the discrepancy between the two results plausibly reflects 
next-order corrections to as well as uncertainties in the determination of A in a neutron 
scattering experiment. 

A different approach is to make use of the scaling result xq ~^^~^) V and neglect 
any possible temperature dependence of a in xq = a^^. To the extent this holds, the full 
temperature dependence of ^(T) may be determined directly from T2G, as shown in Fig.|^. 
Moreover, to the extent that the QNLo" model is valid, despite the absence of a rigorous 



Lorentz invariance of the action, the relation A = c/^ should still approximately apply in 
the QD regime at low temperatures. With these assumptions, for a given choice of a, xq; 
^, and c may be estimated by comparing the neutron scattering result A = 12meV with the 
NMR results for ^^Ti and T2G and assuming Lorentzian form of x{q)- Fo^' convenience, we 
introduce a dimensionless quantity otp = «/«tp7 where otp = 14.8states/eV is a value of 



a estimated by Thelen and Pines [|18| for YBa2Cu306.9. We find that at low temperatures, 
in the QD regime, ^ ~ 4.2(5^? and c ~ a^p ■ 0.19eVA, while in the QC regime, one has 

= (Stp • 0.1(1+T/100K). We note in passing that T-linear QC behavior of both and 
T2G in a doped antiferromagnet has been recently reported by Glenister, Singh, and one of 
us (A.S.) Il; see also Ref. 0. 

We examine next the extent to which the neutron scattering experimental results on 
Lai 85Sro.i5Cu04 are consistent with the NMR results, the magnetic phase diagram of SP, 
and our proposed variation of A with hole concentration. Here there are two issues. First, 
if the hole doping in Lai.85Sro.i5Cu04 is comparable to that found in YBa2Cu306.55, one 
would expect from Fig.^ that the spin gap A for Lai.85Sro.i5Cu04 is ~ 6meV. While the 
experiments reported in Ref. En] show no clear evidence for such a gap, we find an upturn 



near 00 < 6meV in Xl(^) 35K measured by Thurston et al. pT|, as might be expected if 



A ~ 6meV for this material. Second, while the incommensurate spin fluctuation peaks, at 



(tt, tt ± 5) and (vr ± 5, tt), where 5 = 0.2457r, reported in Ref. |2^, are not incompatible with 
our picture, such a large value of 5 is incompatible with NMR experiments. The results of 
Ref. confirm the prediction p3[ that the hyperfine constants for the 214 systems are 



quite close to those inferred for the 123 systems and specifically, for YBa2Cu3 06.63 p4 
On taking the neutron results for x"(q, cu) it is straightforward to show that at lOOK, the 
copper-to-oxygen relaxation rate ratio, (^^Ti/^^Ti y) ^ 35 while the ^^Cu relaxation rate 
anisotropy ratio, (^'^Ti y/^^Ti _|_) ~4.5. These results, which reflect the leakage of correlations 
away from (tt, tt), brought about by the comparatively large value of 5, are to be compared 
to those found experimentally: (1^Ti_||/63Ti_||) >80 ^ and (^^t^ ||/63Ti^^) ~2.6. 

Can the NMR and neutron scattering results be reconciled? The problem is not with the 



magnitude of the measured spin fluctuation peaks. The value of Ti mT is ~ 5.5 x 10 K-sec 



at T = 35K, when calculated from the neutron scattering results of Ref. for q-integrated 
intensity xl(i^), assuming commensurability. This compares favorably with the NMR result, 
^^Tij|T~5.9 • 10~^K ■ sec; indeed the slightly larger value at 35K may reflect the continued 
upturn associated with a spin gap ~ 6meV. The problem rather lies in the assumption that 



the peaks are incommensurate. If instead, as suggested independently by Slichter |]26| and 
Phillips what is being observed is discommensuration, associated with the formation of 
domains in the Cu-0 plane, the NMR and neutron scattering experiments can be reconciled, 
as shown in recent calculations by Monthoux Thus in neutron scattering one is likely 
seeing both commensurate spin fluctuation peaks and domain structure. Two other experi- 
mental results lend support to this proposal: the appearance of a slightly incommensurate 
peak in YBa2Cu306.6 [[H] which might be expected to show behavior not very different from 
that of Lai.85Sro.i5Cu04, and the appearance of a second ^^Cu resonance line in both the 
Sr-doped 214 systems and La2CuO4.032 ||3^, which would plausibly be associated with 



^^Cu nuclei located in or near domain walls. 

In summary, we have shown that NMR, bulk susceptibility, and neutron scattering data 
in the underdoped materials are consistent with the picture of Sokol and Pines where 
the spin gap onset is associated with the crossover from the QC, z = 1, regime to the 
low temperature QD regime. The spin gap magnitudes in YBa2Cu306.63, YBa2Cu408, 
and Lai,85Sro.i5Cu04 compounds determined from different experiments are consistent with 
each other. The magnitude of the spin gap gradually increases as doping increases; since 
the quasiparticle contribution also increases, this eventually leads to the smearing of the 
gap and a crossover to the overdamped, z = 2, regime in the fully doped materials. Such a 
doping dependence is consistent with the theoretical prediction that A~c/.^, provided the 
correlation length decreases more rapidly then c as the doping increases. 
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FIGURES 



FIG. 1. The spin susceptibility for: □ - La2Cu04 |31], from the bulk susceptibility; ■ - 
Lai.85Sro.i5Cu04 ||], from Cu Knight shift; • - YBa2Cu308 ||, from Cu Knight shift; O {^^Kab), 
A C''^Kax), V i^'^Kc), O C''^Kiso) - YBa2Cu306.63 0, from the Knight shift. Solid lines show linear 



fits to the high temperature parts of the respective data, except for YBa2Cu306.637 where the slope 
could not be determined because little high temperature data is available; it has been chosen to be 
the same as for YBa2Cu408 instead. 

FIG. 2. The values of the spin gap in YBa2Cu30j, at different oxygen concentrations y deter- 
mined in inelastic neutron scattering experiments, see the text for discussion; ■ - M. Sato et al. 
|ll, O - J. Rossat-Mignod ei a/ fill, A - J. M. Tranquada a/ [O, O - H. A. Mook et al M\. 



FIG. 3. The Gaussian spin-echo decay time for YBa2Cu306.63) after Ref. Q. The solid line is 
the high temperature linear fit (QC regime); the dotted line represents its saturation in the QD 
regime. An alternate axis shows the temperature dependence of as given by the scaling relation 

r2G~i/^. 
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